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Abstract 

A huge family of solvable potentials can be generated by systemat- 
ically exploiting the factorization (Darboux) method. Staarting from 
the free large class of the known solvable families is thus repro- 

duced, together with new ones. We explicitly find and solve several 
new singular potentials obtained by iteration from the V = case; 
some of them have an E = bound state and constant phase shift 
without being explicitly scale invariant. The new potentials are ra- 
tional functions, and can be related to rational solutions of the KdV 
family. 
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1 Introduction 



The search for solvable potentials in quantum mechanics is an old and large 
industry. In this paper we use the factorization method of Darboux ([1],[2]) 
to formally construct infinite families of fully solvable potentials, all related 
to the free V(x) = case. 

The essence of the method is the following. We start from a hamiltonian 
H and a particular solution (E , <j> ) (D = will be used troughout the 
paper) 

Ho = -D 2 + V (x); HoMx) = EoMx) ( L1 ) 

The solution O needs not to be physical, i.e. it might blow up at finite 
or infinite distances. Then we construct the partner potential 

V 1 (x)-E = W' 2 (x) + W"(x), with 4> {x) = exp(-W(x)) (1.2) 
such that the new hamiltonian H 1 = —D 2 + Vi(x) has as solutions 

H^k = E k ^p k where ^ k = A(p k , A = D + W'(x) and H (p k = E k <p k (1.3) 

In other words, all the solutions (<f) k ,E k ) of the first problem generate 
solutions (A(j>k, E k ) of the new problem. But as A(f) = 0, the new solution is 
ipo = 0o 1 ■ At times, some solutions have to be excluded by physical reasons. 

We shall need also the second solution from the known one ipi at the same 
energy; it is 

i>2 = i>i J ^ 2 dx (1.4) 

These resuls are fairly well-known; for convenience of the reader we supply 
simple proofs of the above statements in Appendix I. 

Our program is to start with the zero potential (free particle) V(x) = 
Vo = and iterate new potentials from its solutions. There are four types of 
different potentials, from four solutions as follows: [3] 

(E = k 2 > 0, k = 1) : = cos(x),^ Vl(x) = +2sec 2 {x) 

(£7 = 0): =x,^V 1 (x) = \ 

x A 
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(E = -k 2 < 0, k = 1) : = cosh(:r),^ Vi(x) = -2cosh~ 2 (a;) 

= sinh(x), =>- Vi(x) = +2 sinh~ 2 (a;)(1.5) 



Notice the last three "solutions" are unphysical; a fifth potential V\ (x) = 
2csc 2 (x) is just the first one displaced |. 

The complete solutions for these first-step potentials are obtained by 
pulling from the free solution by the corresponding A operator, see (1.3). 
Here we recall only the situation for one of the more important cases, provid- 
ing the only case with a potential valid in the whole straight line (— oo, +oo) 

V(x) = -2cosh- 2 (x), W'(x) = - tanh(x) (1.6) 
There is a ground state 

4>o — 0o 1 ( x ) = cosh _1 (x) (unnormalized) (1.7) 
and scattering solutions 

^ki.^) = (D — tanh(x)) exp(ikx) = {ik — tanh(x)) exp(ikx) (1.8) 

corresponding to a transparent (reflectionless) potential, with pure transmis- 
sion 

/ 7 \ ^ / \ 

t(k) — transmission = — (1.9) 

v ' ik + l v ' 

This potential is critical, having an E = resonance; in fact, all trans- 
parent potentials are critical [4] ; our potential (1.6) corresponds to a single 
soliton. 



2 Study of a second-step potential 

The power of the method can be seen now; as the V(x) = case is triv- 
ially solvable for any energy E, physical or unphysical, we have now four 
potentials, all fully solvable; and from each solution of each energy of each 
potential, we can in principle obtain a new, still fully solvable potential. In 
this paper we shall elaborate only in the families associated to the E = 
(intermediate) case. 

We start now from the "centrifugal" potential V(x) = \ and recall the 
(unphysical!) solutions: one is ^ , as <f>o(x) = x is the starting solution for the 
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Figure 1: The potential (2.10) , with two independent parts (// = 1) 



V(x) = 0, E = case; the other is (cf. (1.4)) \ J x 2 dx ~ x 2 . Although both 
solutions blow up at x = and x = oo respectively, they are instrumental in 
obtaining a one-parameter family of bona-fide , physical potentials: 

From the general E = wavefunction (j)(x) = ^+bx 2 , with W'{x) = 
and \i = | > we get the new, second-step potential family 

V,(x) = W'\x) + W'\x) = 6 f x f~^ ] (2.10) 

The new potential(s) is singular: it has a double pole at x — — c , c = 
/i 1 / 3 , < /i < oo. For a; > — c it has an atractive part, and a repulsive tail; 
for x < — c is purely repulsive; both iterpolate between 

V(x) ~ 2 (x ~ -c) ... V(ar) ~ — ^ (x -»• ±oo) (2.11) 
[x -\- C) x 

See Fig.l 

Both V\ — \ and our potential (2.10) correspond to some rational solu- 
tions of the KdV equation [5] ; the relation is interesting and we elaborate on 
it in Appendix II. 
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Of course, the barrier at x = — c is impenetrable: we have two different 
physical problems. 

1) Case x > — c. A bound state candidate with E = —k 2 < would 
behave like exp(— kx) at large x, so we try the E < (unphysical) solution 
4>{x) = exp(-nx) of the V(x) = case and prolongate it twice , as explained 
above. We find 

^oW=M^) = (D- 2 f~/ )(D--)exp(-«x) 

X{X 6 + /jL) X 

= (k 2 + 3x ^p^r) exp(-ra) (2.12) 

For this wavefunction to be physical it has to be zero at the singularity: 
this leads to the eigenvlue equation 1 + k(x = — c) = 0. Hence there is a 
single bound state with energy E = —k 2 , k = 1/c = /i -1 / 3 , and whose 
(un-)normalized wavefunction is 

M*) = (l/c 2 + 2 3X/ * 2 ) exp(-x/c) = { * ,C+l }\ exp(-x/c) (2.13) 
x z — cx + — cx + & 

which behaves in the expected way for a ground state: nodeless, normalizable, 
decaying fast at x — > oo. By construction, this state is t/ie only bound state. 
For E = k 2 > we have total reflection; we write the wavefunction as 

4>k(x) = A 2 Ai(aexp(ikx) + bexp(—ikx)) (2-14) 

and impose ipk(x = — c) = 0; this fixes a/b as a phase, 

a/6 = ~ 1 + ikc exp(2ikc) (2.15) 

From the asymptotic behaviour we extract the S-matrix as usual in scat- 
tering in one radial dimension 

ip k (x » 0) = exp(-ik(x + c)) - S(k) exp(ik(x + c)) (2.16) 

and comparing with (2.14), we derive 

S(k) = \^ (2.17) 
v ' 1 + ike 



5 



Or, for the phase shift S(k) = exp(2i8(k)) 

5(h) = — arctan(/cc) mod n (2-18) 

wich has to be interpreted carefully: with centrifugal tails V(x) — > ^ , /or x » 
the usual rule 5(oo) = does not apply. The interpretation of (2.18) is as 
follows: 

At k — , the bound state contributes +tt to the phase shift (Levin- 
son's theorem) and the long tail (x >> 0) of the potential, which is 6/x 2 = 
l(l + l)/x 2 (1 = 2), contributes -2(tt/2); hence, S(k = 0) = . At very 
large k, the phase shift is dominated only by the short tail, still centrifugal 
+2/(x + c) 2 , which should produce a —tt/2 shift. All this is reproduced by 
(2.18) with the determination arctan(O) = 1 . 

Notice the S-matrix (2.17) is about the simplest with the pole at the 
bound state k = +i/c : this is very similar to the forward amplitude for 
the solitonic scattering (1.9): it seems that the fact that there is a single 
bound state determines the phase shift, and other features of the potential 
are somehow irrelevant. 

2) Case x < — c. Here there is also total reflection, but obviously no 
bound state, and an analogous calculation gives the S-matrix as inverse of 
the previous one, and we get 

S(k) = + arctan(/cc) mod n (2-19) 

At k — the long tail contributes —2 (tt/2), hence we determine arctan(O) = 
— 7r ; as k — > oo , the short tail dominates with S(oo) = —tt/2 ; of course, the 
only invariant statement is the difference, that is, the span A = 5(0) — 5(oo). 

A surprising property of the potential (2.10) has to do with the golden 
ratio $ = (1 + v^5) / 2 : for x > , the maximum xm and minimum x rn of 
V(x) in (2.10) are 

xl = (2 + 3$), x 3 M = (2- 3/$) -1/$) (2.20) 

and the same happens for the values of the potential: 

V{Xm) = 2 ^\l+% 3 ' V{Xm) = V{Xm) (0 (2 ' 21) 

1 It is well known, e.g. in 3D scattering, that a purely centrifugal potential V cen t — 
1(1 + l)/x 2 produces a negative constant phase shift 5 cen t(k) = —lir/2 
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Figure 2: The three regions of the potential (3.22) with (// = 1) 

While we do not fully understand this relation, we notice the same thing 
appears in the KdV for two solitons with velocities k\ and h?: $ = fo/fci, sep- 
arates the overlapping and non-overlapping profiles ([5], p. 190; the discovery 
seems due to Lax) ; it is another intriguing connection between solitons and 
special potentials. 



3 Some generalizations 

For the next step we start with the potential V(x) = 6/x 2 , take the general 
E — solution (f>(x) = fi/x 2 +x 3 and construct, as before, the new, interesting 
potential 

, . 2 6x 10 - 18ux 5 + a 2 . . 

that we plot in Fig. 2 

This potential contains three disconected pieces: 

I- x > . Atraction plus repulsion. 

II c < x < . A confining potential (c = /i 1 / 5 > 0) 

III- x < — c . A repulsive potential. 
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The solutions are again straightforward but tedious, the procedure to 
obtain them is as in the previous case, and we just indicate and quote the 
results: 

l-x > . There is a single bound state with k — and (unnormalized) 
wavefunction 

2 

M*) = = 0" 1 (3.23) 

\i + x b 

and total reflection with wavefunction 

Mx) = A 3 A 2 A l( f) k (x) (3-24) 

with 

A 3 = D + W = D + 2 f " 3X " (3.25) 

x{n + x°) 

and A2 = D — 2/x , A\ = D — 1/x . Hence, if S(k) is the S-matrix for the 
previous V(x) = 6/x 2 and S(k) the new one, 

i; k (x»0) = (D + W'(oo))$ k (x » 0) = (D + W\oo))(exp(-ikx) - S (k) exp(ikx)) 
= iV(exp(— ikx) — S(k) exp(ikx)), hence (3.26) 
S(k) = S{k){W {00) +ik)/{W'{oc)-ik) = {+l){-l) = -l 

because W'(oo) = and S(k) , due to 2(2 + l)/x 2 , is = +1. So 

S{k) = 1 or 5{k) = I mod tt! (3.27) 

These results are worth commenting: First, the E = bound state is 
obvious, because A(f> = in the previous V = 6/x 2 potential implies A^cfi^ 1 = 
0, and _1 zero-less, normalizable. Because of the repulsive tail, it is a 
bona fide bound state, not an E = resonance, so it will contribute +7r to 
Levinson's theorem [7]. 

The constant phase shift is suspitious of some kind of scale invariance. 
In fact, an scale-invariant bound state can exist if at all, at E — ; this is 
our case! The interpretation of the phase shift "span" is this: for k — > 0, the 
bound state contributes +n , the long tail 12/x 2 gives — 37r/2 : so 5(0) = 
— 7r/2 , or S(k = 0) = — 1 . At k — > 00 , the short tail contributes — n/2 , so 
5(k = 00) = — 7r/2 , and the total span of S(k) is zero (while it was +7r/2 in 
the previous case). 
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It is remarkable that a variable (i.e. , not purely centrifugal) potential, 
indeed supporting a (E = 0!) bound state is still "conformal" and produces 
constant phase shift. We offer the following explanation: 

The previous potential V (x) = 6/x 2 is manifestly scale invariant: 

[D, H ] = -2H where H = —D 2 + V = A^A (3.28) 
and D = x ■ D is a dilatation generator. Now 

if i — AA^ — A ■ (A f • A) ■ A' 1 = A ■ H ■ A' 1 (3.29) 

Hence 

[Da, Hx] = -2E X with D A = A- D ■ A' 1 (3.30) 

(Notice A is invertible outside the bound state). Now for x >> , W'(oo) = 
, so A = D + W — > D , and therefore 

Da — > D(x ■ D)D _1 — x-D + 1 — D + l (3.31) 

That is: the traslated symmetry of the new hamiltonian still guarantees 
constancy of the phase shift. 

To the best of our knowledge, this is a first case of a potential, not purely 
centrifugal, with constant phase shift. 

II — c < x < . This confining potential produces of course an uninter- 
esting, infinite ladder of bound states, reminiscent of the potential V(x) = 
— 2cosh _2 (a;) alluded to in Sect. 1. The eigenvalues are E n = k 2 where 

tan (k n c) = — - 3.32 

3 - (k n c) 2 

which is a simple trascendental equation with infinite roots < k\ < k 2 ... < 
k n < ... which tend to nir for n » 1 ; hence the spectrum is asymptotically 
parabolic, as for a particle in an infinite box; this is to be expected, as the 
potential (also in the V(x) = 2sec 2 (x) case) is negligible for higher excited 
wavefunctions. In fact, the normalizable wavefunctions can be written easily, 
but we refrain of doing it. 

III- x < — c . At the left, a purely repulsive potential produces only total 
reflection, and the S-matrix and the phase shift are computed to be 
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So the total span 6(0) — 5(oo) is now = n , and the phase shift is not constant, 
going smoothly from — 3tt/2 to — 7r/2 in the interval k = — > k = oo. Of 
course, "conformal" invariance has been lost because the singular point is at 
x = — c , not at x — 0. 

From the many possible generalizations, we consider in this paper just 
one more case: the general partner of the n-step manifest scale invariant 
potential V (x) = n(n + l)/x 2 . The two E = solutions (both unphysical 
again) are x n+1 and x~ n ; so defining 

<j>(x) = fi/x n + x n+l n > (3.34) 

the partner family is 

(n + l)(n + 2)x 4n+2 - Qfinjn + l)x 2n+1 + /i 2 n(n - 1) 

which again exhibits the three regions as before. In particular 

x > : a partly attractive potential, which supports again just a bound 
state at zero energy ; total reflection occurs with (again) constant phase shift. 
The bound state is <p~ l , of course, and it turns out that 

S(k) = (-1)" +1 (3.36) 

by the same argument as before, namely S(k) = —S n (k) , where S n (k) is the 
S-matrix for Vq(x) = n(n + l)/a; 2 , namely S n (k) = (— l) n . 

We have therefore found and infinite family of "scale" invariant potentials, 
with a unique E — normalizable bound state, and constant (in fact ±1) 
S-matrix. The (modified) Levinson theorem applies; namely the span 5(0) — 
S(oo) is zero: at low k, there is a +n contribution from the bound state, and 
—n(n + l)7r/2 value from the long tail. At large k , the short tail takes over, 
contributing — (n — 1)tt/2 . The constancy of 5(k) comes, as before, from the 
appropiate conjugation of the manifest dilatation symmetry of the previous 
potential, just as in the worked-out case n = 2. 

In the confining region — c < x < , with c = +/j l 1 /( 2n+1 ) > , there is a 
pure point spectrum, with again a limiting parabolic growth in the energy. 
The spectral equation is a natural generalization of (3.32); we state only the 
next case, n = 3; the trascendental eigenvalue equation is 

. 7(kc) 3 - 105(fcc) 
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Finally, in the pure repulsive part of the potential, x < — c , there is only 
total reflection with a simply variable phase shift. The total span is 



5{0) - 6(00) = -{n + l)7r/2 - (-7T/2) = -mr/2 (3.38) 

because the potential behaves like +2/(x + c) 2 close to the pole. The exact 
S-matrix can be calculated as before. We just quote the result only again for 
n = 3 : 

= 42(fcc) 2 + 105-*(7(fcc) 3 -105(A:c)) 
(complex conjugate) 

The general S(k) starts at 5(0) = +1 for n — 3, 5, 7, ... , and 5(0) = — 1 
for n even; it becomes S(oo) = —1 after n half-turns. The phase shift 
connects smoothly — (n + l)n/2 at k = with — tt/2 at k = 00 . 

We can see also why the first case n — 1 is special: at right the potential 
is +6/x 2 for n = 2 , and at x = is V = , as n — 2 = so in this case there 
are only two regions with no confining part. 



4 Other potentials 

Once the general procedure is understood, is a matter of mechanical calcula- 
tions to find and to solve any other V = 0-related potentials. We shall report 
on a full investigation elsewhere [11] 

Here we just report that we can, by our procedure, recover many of the 
"shape invariant" potentials in the review Infeld-Hull paper [1] ; in fact, all 
the families included in the "A-type" classification of [1] . The other types 
B...I are in some way degenerate: they include, among others, the oscillators, 
Kepler and Morse potentials, which are not directly connected to the V = 
case, but still are "shape invariant" and solvable. As shown in [11] , the 
Kepler problem is related to the V = potential in a constant curvature 
(spherical for bound states) space. 

The natural minimal generalization of V = 2/x 2 is obviously the centrifu- 
gal potential 

= n(n + l) n = 012 ___ (440) 

x A 

This is obtained from V = by making use of the solutions Ob j Ob j Ob j • • • j Ob 
in each step. 
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The minimal natural extension of V = 2 sec 2 (x) is 



V(x) = +n{n + l)sec 2 (V) n = 0,l,2... (4.41) 

The intertwining superpotential satisfies W n {x) = ntan(x) , with <f> n (x) = 
cos™ (re) as the generating wavef unction.; notice the energy scale gets dis- 
placed; this confining-potential family contains a pure discrete spectrum, 
approaching the parabolic infinite-box situation. 

Similarly 

V(x) = -n(n + I) cosh- 2 (x) n = 0,l,2... (4.42) 

comes from W^(x) = — ntanh(rr) and (f> n (x) = cosh ra (x). 

There are n bound states and a E = resonance, plus perfect transmis- 
sion (no reflection); it is the well known "n-solitonic" potential, with all the 
elementary solitons on top of each other at x — [5] . 

The final minimal family is 

V(x) = +n(n + 1) sinh- 2 (a;) n = 0,1,2,... (4.43) 

This comes from W^(x) = +ncoth(a;) and (f> n (x) = sinh n (a;) . It corresponds 
to total reflection, with variable phase shifts, and no bound states. 

All these four families are still exactly solved even for n — > A noninteger, 
by "prolongation" (see [1] or the review [8]) ; they are shape-invariant [12] 
and therefore included in [1] . As they are not related directly with the 
vacuum V = case, we do not discuss them. 

The only potential of "A" type of [1] not include so far is 

T _ . a + bcos(x) fAAA ^ 
V(x) = . 2 , y 4.44 

This can still be also obtained in our scheme in the following, indirect way: 
the potential 

V(x) = -V* (4.45) 
sin [x) 

is a prolongation of the Vi(x) = 2 csc 2 (x) of § 1, and it admits the unphysical 
eigenfunction 

(f)(x) = y/sm(x) cot(x/2) (4.46) 
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Hence the corresponding partner potential is, with W'(x) = — cot(x/2) + 
csc(a;) , 

which is of type (4.44). The energy of the unphysical solution is +1/4 , 
whereas the ground state of (4.47) is sin 3 / 2 (x) , with energy = +9/4 . We 
conclude that the "A" type family of Infeld-Hull [1] can be included also in 
our scheme of things. 



5 Conclusion 

The whole set of analytically soluble potentials (not to speak of the quasi- 
soluble ones [13] ) is very, very large. In this paper we have shown how 
starting with the free case, V(x) = , and just by playing around with the 
unphysical solutions for E — only, a large family is obtained; the generic 
case includes a confining potential defined in a segment of the line, a purely 
repulsive half-line defined potential, and an also half-line defined potential, 
supporting a bona fide unique E = bound state with trivial (i.e. constant 
= ±1) S-matrix. 

The natural generalization of the four different potentials obtained in the 
first step from V(x) = includes all the non-degenerate cases in the Infeld- 
Hull series, if we include prolongations, that is, substituting n(n+ 1) n e N, 
by A(A + 1) for arbitrary, real positive A. They correspond to solutions of 
the hypergeometric equation, which is also related to the SL(2, R) group. ; 
the degenerate I-H cases B...I (i.e. Coulomb,...) correspond to solutions of 
the confluent hypergeometric equation. 

There is still work in progress; we have not exhausted even the E = 
family (for example, we can iterate the potential (3.35)!). As stated, we plan 
to report on other cases in a later publication; see also [11] 

The two Appendices explain the Darboux method and elaborate on the 
KdV connection, as promised. 
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A Appendix 

We probe the results used in §1. From = (-D 2 + V)^> = E^f, define 

exp(-W(a;)) = * (A.l) 

W satisfies a Riccati first order equation V — E = W' 2 — W": the scale 
invariance ip — > \ip becomes translation invariance for W, hence W itself 
does not appear in the new equation, ip needs not to be physical, i.e. it can 
be singular. But now the hamiltonian factorizes: 

-D2 + W' 2 (x) - W"{x) = (-D + W\x)){D + W\x)) = A ] A (A.2) 

with A = D + W'(x). We obtain 

A^Am = (H-E)m (A.3) 

The partner hamiltonian is defined as H' = AA^ + E, so 

H<j>' = E'cj)' = (A^A + E)<j>' {AA^ + E')A<j>' = E' A<p' (A.4) 

For each solution (0', E') of the former H we obtain a solution (A<f>', E') 
of the new H' . Of course, 0' might be physically unacceptable. This is the 
essence of the method. 

Now for the second solution. = exp(— W) implies A<p = 0, A = D + 
W. Or (exp(-W) • D ■ exp(+W)) exp(-w) = 0, hence (exp(+W) ■ D ■ 
exp(— W))^ 1 = 0, A^ifi^ 1 = so corresponds to for E' = E. Now 
if 0' is the second solution of the original H with energy E, A^A<p' = but 
Acj)' ^ 0, hence Acp' is in the kernel of A\ and therefore 

exp(-W^) • D ■ exp(+l^)0 / = exp(W), or 0' = J 0" 2 dx (A.5) 

as stated. 

The use of the second solution to generate new potentials seems to start 
with [9] ; see also the previous work of Abraham and Moses [10] 
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B Appendix 



The KdV equation (u = u(x,t), U jt = fjr etc.) 

u :t = 6uu tX - u >xxx (B.l) 

is one of the deformation equations associated to the Schrodinger equation, 
H(fj,) = —D 2 +u(x,fj), where i) is the deformation parameter. For 
this reason some simple solutions of KdV are interesting potentials for the 
linear problem; we take some results from [5] . The travelling wave solution 
u = u(x — vt) 

u(x, t) = -(l/2)v cosh" 2 (v/i;(a; - vt - x )) (B.2) 

corresponds to our first step with E < and = cosh (a;); it is the solitonic 
potential. Multisolitonic potentials correspond to iteration from this solution, 
but these are not considered in this paper. 

Rational solutions of KdV are closer to our potentials; for example u = 
2/x 2 arises as the simplest t-independent rational solution, and it is our first 
potential from the E = 0, = x solution. The natural scaling invariance of 
KdV x — > \x, t — > A 3 t, u — > \~ 2 u leads at once to the rational solution 

, x 6x(x 3 - 24*) 
u ( x t) = —\ -£. B.3 

which is our potential (2.10) (with \x = 12t). Similarly the other rational 
potentials we obtain are connected with rational solutions of higher-order 
KdV-hierarchy equations; we shall report on a full investigation elsewhere. 
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